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We study the effect of isospin degree of freedom on the balance energy (Ebai) as well as its 
mass dependence throughout the mass range for two different sets of isobaric systems with 
N/Z = 1 and 1.4 at different colliding geometries ranging from the central to peripheral 
ones. Our findings reveal the dominance of Coulomb repulsion in isospin effects on Ebai as 
well as its mass dependence throughout the range of the colliding geometry. Our results 
also indicate that the effect of symmetry energy on the energy of vanishing flow is uniform 
throughout the mass range and throughout the colliding geometry. 
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1 Introduction 



The existing and upcoming radioactive ion beam (RIB) facilities [1-4] around the world 
led to the increased interest in isospin degree of freedom and equation of asymmetric 
nuclear matter. The above study helps to isolate the isospin-dependent part of nuclear 
equation of state which is vital in understanding the astrophysical phenomena. Isospin 
degree of freedom plays an important role throiigh both the nuclear matter equation of 
state and as well as via in-medium cross section. Therefore, heavy-ion collisions induced 
by neutron-rich nuclei provide a good opportunity to explore the isospin dependence of 
in-medium nuclear interactions. 

After about three decades of intensive efforts in both nuclear experiments and theo- 
retical calculations, equation of state for isospin symmetric matter is now relatively well 
understood. The effect of isospin degree of freedom on the collective transverse in-plane 
flow as well as on its disappearance [5-7] (there exists a particular incident energy called 
balance energy (Eba/) or energy of vanishing flow (EVP) at which transverse in-plane flow 
disappears) has been reported in the literature [8-10], where it was found that neutron- rich 
systems have higher E^ai compared to neutron-deficient systems at all colliding geometries 
varying from central to peripheral ones. The effect of isospin degree of freedom on Ef,ai 
was found to be much more pronounced at peripheral colliding geometries compared to 
central ones. As reported in the literature, the isospin dependence of collective flow as 
well as its disappearance has been explained as a competition among various reaction 
mechanisms, such as nucleon-nucleon collisions, symmetry energy, surface property of the 
colliding nuclei, and Coulomb force. The relative importance among these mechanisms 
is not yet clear [8]. In recent study, we [11] confronted theoretical calculations (using 
isospin-dependent quantum molecular dynamics (IQMD) model [12]) with the data at all 
colliding geometries and were able to reproduce the data within 5% on the average at all 
coUiding geometries. Motivated by this good agreement of the calculations with data, we 
[11] studied the isospin effects on the Ef,ai throughout the mass range 48-350 for two sets of 
isobaric systems with N/Z — 1.0 and 1.4 at semi central coUiding geometry. These results 
showed that the difference between the E^a; for two isobaric systems is mainly due to the 
Coulomb repulsion. It was also shown that Coulomb repulsion dominates over symmetry 
energy. These flndings also indicated towards the dominance of the Coulomb repulsion in 
larger magnitude of isospin effects in E^a; at peripheral collisions. Here we aim to extend 
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the study over full range of colliding geometry varying from central to peripheral ones. 
Section 2 describes the model in brief. Section 3 explains the results and discussion and 
Sec. 4 summarizes the results. 



2 The model 

The present study is carried out within the framework of the IQMD model [12] which 
is an extension of the QMD model [13-15]. The IQMD model treats different charge 
states of nucleons, deltas, and pions explicitly, as inherited from the Vlasov-Uehling- 
Uhlenbeck (VUU) model. The IQMD model has been used successfully for the analysis 
of a large number of observables from low to relativistic energies. The isospin degree of 
freedom enters into the calculations via symmetry potential, cross sections, and Coulomb 
interaction. 

In this model, baryons are represented by Gaussian-shaped density distributions 
fi{f,p,t) = ^;^cxp(-[f-rl(t)]2^) 

xeM-\p-m?^) (1) 

Nucleons are initialized in a sphere with radius R = 1.12 A^/^ fm, in accordance with 
liquid-drop model. Each nucleon occupies a volume of h^, so that phase space is uni- 
formly filled. The initial momenta are randomly chosen between and Fermi momentum 
(pf)- The nucleons of the target and projectile interact by two- and three-body Skyrme 
forces, Yukawa potential. Coulomb interactions, and momentum-dependent interactions. 
In addition to the use of explicit charge states of all baryons and mesons, a symme- 
try potential between protons and neutrons corresponding to the Bethe-Weizsacker mass 
formula has been included. The hadrons propagate using Hamilton equations of motion: 

dfi d{H) dpi d{H) 



dt dpi ' dt dfi 



(2) 



with 



(H) = {T) + {V) 

i i j>i 

xfj{f',p',t)dfdf'dpdp'. (3) 
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The baryon potential , in the above relation, reads as 



= [t,S{f'-r} + t2S{f'-r}p'^-'C-^)] 

exp(|(f^-f)|//x) Z,Z,e^ 
(|(r'-f)|//i) '■|(r-''-f)| 
+i4ln'[i5(p'-p)' + l](^(r '-r) 

+t,-Ts,Tsj5{fi'-fj). (4) 

f?o 

Here -Zi and Zj denote the charges of ith and jth baryon, and T^i and T^j are their re- 
spective components (i.e., 1/2 for protons and —1/2 for neutrons). The parameters 
II and ^i,....,Z6 are adjusted to the real part of the nucleonic optical potential. For the 
density dependence of the nucleon optical potential, standard Skyrme type parametriza- 
tion is employed. Note that such Skyrme forces and related potential are essential tool to 
understand low energy phenomena such as fusion, fission and cluster radioactivity [16]. 
The momentum-dependence V^^- of the nn interactions, which may optionally be used in 
IQMD, is fitted to the experimental data in the real part of the nucleon optical potential. 

3 Results and discussion 

For the present study, we simulate several thousands events of each reaction at incident 
energies around E^ai in small steps of 10 MeV/nucleon. In particular, we simulate the reac- 
tions 2^Mg+24Mg, 58Cu+5^Cu, 72Kr+^2j^r, ^^Cd+^^Cd, i20Nd+i20Nd, ^^^Eo+^^^Ro, hav- 
ing N/Z = 1.0 and reactions ^^Ne+^^Ne, ^^Ct+^^Ct, ^^Zn+^^Zn, ^^Zr+^^Zr, i^osn+i^ogn, 
and ^^^Ba-|-^^^Ba, having N/Z = 1.4, respectively, in the whole range of colliding geome- 
try. The colliding geometry is divided into four impact parameter bins of 0.15 < 6 < 0.25 
(BIN 1), 0.35 < b <0.45 (BIN 2), 0.55 < 6 < 0.65 (BIN 3), and 0.75 <b< 0.85 (BIN 4), 
where b = h/hmax- 

Here N/Z is changed by keeping the mass fixed. We use anisotropic standard isospin- 
and encrgy-dcpcndcnt nn cross section a = 0.8 a^N^^'^'^- The details about the clastic and 
inelastic cross sections for proton-proton and proton-neutron collisions can be found in 
Ref. [12]. The cross sections for neutron- neutron collisions are assumed to be equal to the 
proton-proton cross sections. We also use soft equation of state along with momentum- 
dependent interactions (MDI). The results with the above choice of equation of state and 
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System Size (A) 

Figure 1: (Color online) Ef,ai as a function of combined mass of system for different impact 
parameter bins. Solid (open) symbols are for systems having N/Z = 1.0 (1.4). The solid 
(dashed) lines are power law fit a A"^. r values for full Coulomb calculations are displayed 
in figure. The values of r for reduced Coulomb calculations are given in the text. 



cross section were in good agreement with the data. The reactions are followed until the 
transverse flow saturates. The saturation time varies form 100 fm/c for lighter masses to 
300 fm/c for heavier masses. For transverse flow, we use the quantity " directed transverse 
momentum {pf '')" which is defined as [17] 

{pT) = jJ2-''9HympA^), (5) 

2=1 

where y{i) is the rapidity and Px{i) is the momentum of i*^ particle. The rapidity is 
defined as 

yW = i,„||!i±Ml, (6) 

where E(i) and p^i''-) respectively, the energy and longitudinal momentum of i*'* 
particle. In this definition, all the rapidity bins are taken into account. A straight line 
interpolation is used to calculate E^ai- 

Figure 1 displays the mass dependence of Ebai for four impact parameter bins. The 
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Figure 2: (Color online) r as a function of reduced impact parameter. Lines are only to 
guide the eye. 

solid (open) green circles indicate E^ai for systems with lower (higher) N/Z. Lines are 
power law fit a A^. Ebai follows a power law behavior cx for both N/Z = 1 and 1.4 
(r being labeled as ti q and ri.4 for systems having N/Z = 1 and 1.4, respectively) at 
all colliding geometries. Also, for the clarity, we did not include the very light system 
(A=48) as isospin effects are negligible at all colliding geometries for this system. 

From fig. 1 we see that isospin effects are clearly visible for all the four bins as neutron- 
rich system has higher Ehai throughout the mass range in agreement with the previous 
studies [8,9, 11]. The magnitude of the isospin effects increases with the increase in the 
mass of the system at all colliding geometries. The effect is much more pronounced at 
larger colliding geometries. One can see that the difference between ri.o and ri.4 increases 
with increase in the impact parameter. In Ref. [18], one of us and coworkers studied the 
isospin effects on the mass dependence of E^ai for BIN 2. There the Coulomb potential was 
reduced by a factor of 100 and was showed that the Coulomb repulsion plays dominant 
role over symmetry energy in isospin effects on Ef,ai as well as its mass dependence at semi 
central colliding geometry (BIN 2). Since here we plan to extend that study over a full 
range of colliding geometry, so here also we reduce the Coulomb potential by a factor of 
100 and calculate the Ebai throughout the mass range at all colliding geometries. Solid 
(open) diamonds represent Ebai calculated with reduced Coulomb for systems with lower 
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Figure 3: (Color online) E^ai as a function of impact parameter for different system masses. 
Symbols have same meaning as in Fig.l. Lines are only to guide the eye. [18] 

(higher) neutron content. Lines are power law fit a A^. The values of ti.q (tia) are -0.17 
(-0.14), -0.17 (-0.19), -0.24 (-0.26), and -0.31 (-0.29) for BIN 1, BIN 2, BIN 3, and BIN 
4, respectively. Interestingly, we find that the magnitude of isospin effects (difference in 
Efyai for a given pair) is now nearly same throughout the mass range which indicates that 
the effect of symmetry energy is uniform throughout the mass range. This is true for all 
the colliding geometries. This is supported by Ref. [17] where Sood and Puri studied 
the average density as a function of mass of the system (throughout the periodic table) 
at incident energies equal to E^ai for each given mass. There they found that although 
both Efyai and average density follows a power law behavior cc A'^ , Ehai decreases more 
sharply with the combined mass of the system (with r = -0.42), whereas the average 
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density (calculated at incident energy equal to Ehai) is almost independent of the mass of 
the system with r = -0.05. It is worth mentioning here that the trend will be different 
at fixed incident energy in which case density increases with increase in the mass of the 
system [19,20]. Prom fig. 1 one can also see that the enhancement in E^a^ (by reducing 
Coulomb) is more in heavier systems as compared to lighter systems for all colliding 
geometries. The effect is more pronounced at higher coUiding geometries. Moreover, 
throughout the mass range at all colliding geometries, the neutron-rich systems have less 
Ehai as compared to neutron-deficient systems when we reduce the Coulomb. This trend 
is quite the opposite to the one which we have when we have full Coulomb. This (as 
explained in Ref. [18] also) is due to the fact that the reduced Coulomb repulsion leads to 
higher E^ai- As a result, the density achieved during the course of the reaction will be more 
due to which the impact of the repulsive symmetry energy will be more in neutron-rich 
systems, which in turn leads to less E},ai for neutron-rich systems. 

In fig. 2, we show the variation in At as a function of colliding geometry for both 
full and reduced Coulomb where At = T1.4 — ti.q. Half filled circles (diamonds) represent 
calculations with full (reduced) Coulomb. At values are plotted at upper limit of impact 
parameter for each bin. Lines are only to guide the eye. Prom figure (green circles), we see 
that the difference between T1.4 and ti.o increases with increase in impact parameter which 
is due to the sharp increase in value of |ri|, whereas when we reduce the Coulomb (orange 
diamonds), this difference is negligible, which shows an enhanced effect of Coulomb at 
higher impact parameters. 

In fig. 3a, 3b, and 3c, we display Ef,ai as a function of b for masses 116, 192, and 240, 
respectively, for both full and reduced Coulomb. Symbols have the same meaning as in fig. 
1. Por full Coulomb (green circles), for all the masses at all colliding geometries, system 
with higher N/Z has larger E^ai in agreement with previous studies [8,9]. Moreover, the 
difference between E^ai for a given mass pair, increases with increase in colliding geometry. 
This is more clearly visible in heavier masses. Also for N/Z = 1.4, E^ai increases with 
increase in impact parameter as expected [21,22]. However for N/Z = 1, this is true only 
for lighter mass system such as A = 116. For heavier masses E^ai infact begins to decrease 
with increase in impact parameter in contrast to the previous studies [8,9,11,21,22]. 
However, when we reduce the Coulomb (by a factor of 100 (diamonds)), we find that: 

(i) Neutron-rich systems have smaller E^ai as compared to neutron-deficient systems as 
mentioned previously also. This is true at all the colliding geometries throughout the mass 
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Figure 4: (Color online) The time evolution of number of collisions for various system 
masses at 50 MeV/nucleon for BIN 1 and BIN 4. Various lines are explained in the text . 

range. This clearly shows the dominance of Coulomb repulsion over symmetry energy in 
isospin effects throughout the mass range at all colliding geometries. 

(ii) The difference between for systems with different N/Z remains almost constant 
as a function of colliding geometry which indicates that the effect of symmetry energy is 
uniform throughout the range of b as well. This also shows that the large differences in 
Ebai values for a given isobaric pair are due to the Coulomb repulsions. 

To see the role of Coulomb repulsion on nn collisions, in fig. 4 we display no. of 
collisions for A = 48, 116, and 240 for calculations with and without Coulomb potential 
at central (BIN 1) and peripheral (BIN 4) colliding geometry at incident energy of 50 
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MeV/nucleon. Solid (dashed) lines represent Coulomb full (reduced) calculations. Higher 
(lower) lines represent the results for BIN 1 (BIN 4). We find that Coulomb decreases the 
no. of collisions in medium and heavier mass systems for central bin whereas for peripheral 
bin, the effect of Coulomb on collisions is only for heavier masses (lower panel). 

4 Summary 

We have studied the isospin effects in the disappearance of flow as well as its mass depen- 
dence throughout the mass range for two sets of isobaric systems with N/Z = 1 and 1.4 
in the whole range of colliding geometry. Our results clearly demonstrate the dominance 
of Coulomb repulsion in isospin effects on Ebai as well as its mass dependence throughout 
the range of coUiding geometry. The above study also shows that the effect of symmetry 
energy on E^ai is uniform throughout the mass range and colliding geometry. 

This work has been supported by a grant from Indo-French Centre For The Promotion 
Of Advanced Research (IFCPAR) under project no. 4104-1. 
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